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NOMENCLATURE 


a.,,b..  Arbitrary  interlaminar  shear  coefficients 

1  j  ij 

a,b  Vectors  of  arbitrary  interlaminar  shear 

coefficients 

b  Plate  width 

E  ,E  Longitudinal  and  transverse  Young's  moduli  of 

£k  k  k'th  lamina,  respectively,  with  respect  to 

material  axes 

e  Subscript  denoting  linear-elastic-material 

component 

e,?,g  Vectors  of  arbitrary  displacement  coefficients 

P'  Stress-energy -density  function 

Elastic,  elasto-plastic,  and  total  stress- 
energy-density  functions  for  the  k’th  lamina  as 
a  function  of  lamina  inplane  stresses, 
respectively 


Elastic,  elasto-plastic,  and  total  stress- 
energy  -density  functions  for  the  matrix  material 
between  k'th  and  (k+l)st  lamina  as  a  function  of 
transverse  stresses  acting  between  those  two 
layers,  respectively 


G  Matrix  material  shear  modulus 

m 

G  k'th  lamina  inplane  shear  modulus  with  respect 

Lk  to  lamina  material  axes 

h  Laminate  thickness 

hp(p=l,2,3)  Vectors  of  arbitrary  stress  coefficients 


f.,  N-dlmensional  unit  vector 

N 

K  ,Kt  ,K  t  Ramberg-Osgood  type  constants  for  k'th  lamina 

\  k  k  relative  to  the  material  axes 

k  Index  representing  the  k'th  lamina 
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Plate  length 


Total  number  of  laminae  in  the  laminate 


N  ,N  ,N 
x*  y  xy 


Laminate  inplane  stress  resultants 


n  ,n^  ,n 

fck  £tk 


Ramberg-Osgood  type  exponents  for  the  k'th 
lamina  relative  to  the  material  axes 


P,  (1=1,2, 3) 
xk 


Laminate  membrane  stress  resultant  vector 

Subscript  denoting  nonlinear-elastic-material 
component 

Exponents  analogous  to  Ramberg-Osgood-type 
exponents  defined  in  Reference  17 

Ramberg-Osgood-type  exponent,  defined  in 
Reference  17 

Heat  flux 

Superscript  denoting  matrix  transpose 
Laminate  thickness 

k'th  lamina  thickness  over  which  the  extensional 
stresses  o  and  a  and  the  shearing  stress 


Thickness  of  matrix  material  between  the  k'th 
and  (k+l)st  lamina's  median-surfaces  over  which 
the  interlaminar  shear  stresses  t  and  T 


Uk’Vk 


Modif ied-Reissner  functional 

Median-surface  x  and  y  displacement  functions  of 
the  k'th  lamina,  respectively 

N-dimenslonal  x  and  y  displacement  vectors, 
respectively 

Volume 


Lateral  displacement  funciton  of  laminate 


x,y,z 


Plate  coordinates 
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(1 

xy  ** 


y 


xy 


k 


Y 


eef  f 


k 


Laminate  averaged  inplane  shear  stresses 
relative  to  plate  axes 

Median-surf ace  shearing  strains  of  k'th  lamina 
relative  to  plate  axes 

Interlaminar  shearing  strains  acting  between  the 
median-surfaces  of  the  k'th  and  (k+l)st  laminae 
relative  to  plate  axes 

Effective  strain  defined  in  Equation  ( 1 8 ) 


n12,2,n12,l 

\ 

nt ,  it  *  n  i.  Jit 

> 

%t,t’ nit,  i 


Lamina  coefficients  of  mutual  influence  relative 
to  material  axes 


Coefficients  of  mutual  Influence  relative  to 
material  axes 


0 


Temperature 


IC 


Thermal  diffusivity 


v 


Major  and  minor  Poisson's  ratios  of  k'th  lamina, 
respectively,  relative  to  the  material  axes 


o 


» 


Median-surface  extensional  stresses  of  k'th 
lamina  relative  to  the  material  axes 


o 


Median-surface  extensional  stresses  of  the  k'th 
lamina  relative  to  the  plate  ares 


\ 

Laminate  averaged  extensional  stresses  relative 
>  to  plate  axes 


Median-surface  shearing  stress  of  k'th  lamina 
relative  to  material  axes 


Median-surface  shearing  stress  of  k'th  lamina 
relative  to  plate  axes 
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(0)  (1) 

Txy  'Txy  •** 


Laminate  averaged  inplane 
to  plate  axes 


shear  stress  relative 


T 


yz 


Interlaminar  shear  stress 


vectors 


r 


(0) 

yz 


1 


•  •  • 


Laminate  averaged  transverse  shearing  stresses 
relative  to  plate  axes 


Anamet  Laboratories  Inc. 

HAYWARD.  CALIFORNIA 

PROJECT  SUMMARY 


CONTRACT  NUMBER:  F33615-86-C-321 8 

CONTRACTOR:  Anamet  Laboratories,  Inc. 

3400  Investment  Blvd. 

Hayward,  CA  94545-3811 

PRINCIPAL  INVESTIGATOR:  Rocky  Richard  Arnold,  Ph.D. 

TITLE  OP  THE  PROJECT:  Rapid  Thermal  Loading  of  Delaminated 

Composite  Structures 


TECHNICAL  DESCRIPTION 

The  research  program  proposed  herein  was  directed  toward  develop¬ 
ing  a  theoretical  approach  to  the  problem  of  rapid  thermal  loading 
of  delaminated  composite  structure.  The  structures  examined 
include  flat  and  shallow-curved  plates  and  cylindrical  shells. 
Using  Hamilton's  principle  and  the  Reissner  variational  theorem, 
a  new  dynamic  thermoelastic  variational  principle  was  developed. 
Application  of  this  principle  to  both  plates  and  shells  provided 
the  equations  of  dynamic  equilibrium  and  the  associated  natural 
and  geometric  boundary  conditions.  Important  physically  observ¬ 
able  characteristics  of  composites,  not  usually  found  in  clas¬ 
sical  approaches,  were  also  included  in  the  theory.  These 
characteristics  include  transverse  shear,  material  nonlinearity 
and  delamination  mechanisms. 

In  this  Phase  I  effort,  the  therraoelastlc  model  was  defined 
and  the  governing  equations  derived  for  both  contiguous  and  de¬ 
laminated  plate  and  shell  structures.  Solution  of  the  governing 
equations  and  numerical  examples  constitutes  the  major  emphasis 
of  the  Phase  II  work. 


COMMERCIAL  APPLICATION 

The  most  Immediate  benefit  from  the  proposed  effort  will  be  the 
ability  to  design/analyze  delaminated  composite  structures  with 
improved  resistance  to  rapid  therm3 1  loading.  This  his  direct 
application  to  various  military  systems,  especially  aircraft  that 
may  be  subject  to  hostile  laser  threats.  Other  potential  appli¬ 
cations  include  future  hypersonic  vehicles  and  civilian  transport 
aircraft,  wherein  lower  rates  of  thermal  loading  are  of  concern. 
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1.0  INTRODUCTION 

This  Is  the  Phase  I  final  report  for  work  accomplished  under 
U.S.  Air  Force  Contract  No.  F33615-86-C-3218,  as  a  part  of  the 
Federally-sponsored  Small  Business  Innovation  Research  (SBIR) 
program.  The  period  of  technical  performance  was  from  July  31, 
1986  to  January  31,  1987.  The  research  and  development  effort 
documented  herein  was  accomplished  for  AFWAL/GLXPF. 

The  primary  objective  of  the  Phase  I  research  was  to  derive 
a  set  of  governing  equations  for  both  contiguous  and  delaminated 
composite  plates  and  shells  subject  to  intense  rapid  thermal 
heating.  These  equations  are  amenable  to  programming  on  a  digital 
computer,  which  is  the  primary  emphasis  of  any  potential  Phase  II 
work.  The  objective  of  this  research  has  been  successfully 
accomplished  and  the  assumptions  and  derivations  used  to  complete 
this  project  are  documented  within  this  final  report. 


2.0  THKORY 

The  basic  structural  configurations  under  consideration 
herein  a-e  the  flat  plate,  shallow-curved  plate,  and  complete 
circular  shell.  Schematic  diagrams  of  these  basic  aircraft 
structures  and  the  coordinate  systems  used  in  the  ensuing  deri¬ 
vations  are  included  in  Figure  1.  For  all  structural  configura¬ 
tions,  the  laminated  composite  is  examined  in  two  distinctly 
different  states;  that  is,  either  the  composite  is  undamaged 
(contiguous  state),  or  the  composite  is  delaminated  presumably 
through  the  action  of  a  rapid  thermal  pulse.  In  both  cases,  the 
primary  theoretical  derivations  are  identical--dif ferences  between 
either  the  contiguous  or  delaminated  composite  configurations  are 
manifested  by  the  appropriate  choice  of  stress  and  displacement 
functions  which  will  be  explained  in  more  detail  in  a  latter  part 
of  this  report  (Section  3«1)* 

In  this  section,  the  basic  theoretical  derivations,  valid 
for  both  contiguous  and  delaminated  composite  configurations,  are 
developed  and  limiting  assumptions  stated. 


1 


Figure  1  Curved  Element  Geometry  and 


Anamet  Laboratories,  Inc. 

HAYWARD,  CALIFORNIA 


2.1  Statement  of  Problem  and  Principal  Assumptions 

The  plate  and  shell  elements  consist  of  N  homogeneous  and 

materially  nonlinear  anisotropic  laminae  with  each  lamina  having 

thickness  t„  .  Each  element  of  a  lamina  midsurface  can  undergo 
fmk 

three  translations ,  two  inplane  displacements,  u  and 

1C 

v  (k-1,2, . . .N),  and  a  lateral  displacement  w  which  is  common  to 

all  laminae.  The  inplane  stresses  a  ,  a  and  T  are  assumed 

xk  yk  xyk 

constant  through  the  thickness  of  each  lamina. 

The  bending  stiffness  of  an  individual  lamina  is  assumed 

negligible  compared  to  that  of  the  laminate;  that  is,  laminae  are 

assumed  to  behave  as  membranes.  Consequently,  transverse  shear 

effects  are  accommodated  by  permitting  relative  movement  between 

the  median  planes  of  adjacent  laminae.  The  matrix  material 

between  adjacent  laminae  midsurfaces  is  assumed  to  carry  all  of 

the  transverse  shear  in  the  y-z  and  z-x  planes.  The  normal 

stress  (o  )  is  taken  to  be  negligible  in  comparison  with  the 
z 

other  interlaminar  stresses  (x  ,x  )  for  the  composites  plates 

yz  zx 

considered  herein;  namely,  those  that  buckle  in  the  linear-elastic 

range  and  become  elasto-plastic  in  the  postbuckling  range. 

The  assumptions  cited  above  are  known  to  be  valid  for  the 

contiguous  composite  [1];  however,  for  a  delaminated  composite 

the  normal  stress  (0  )  that  would  exist  at  the  interface  between 

z 

adjacent  laminae  that  contain  the  delamination  is  not  necessarily 
negligible.  During  the  Phase  II  work,  a  closer  examination  of 
this  assumption  will  be  accomplished  and,  if  required,  suitable 
modifications  made  to  the  approach  proposed  herein. 
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2.2  Strain-Displacement  Relations 

The  ultimate  success  of  any  theoretical  model  is  the 
appropriate  definition  of  initial  assumptions — these  initial 
assumptions  determine  the  accuracy  and  utility  of  the  finally- 
developed  analytical  method.  For  this  research,  the  membrane 
strain-displacement  relations  of  Sanders  [2]  are  used.  For  a 
circular  cylindrical  shell,  these  relations  are: 


e  =  —  +  -if— 
x  ax  2laxJ 


(l) 


av  w  l  r  aw  v  ,2 
ey  ay  "  R  2  ^ay  R ' 


(2) 


av  au  aw  aw  +  v  aw 
Yxy  ax  +  ay  +  ax  ay  +  R  ax 


(3) 


For  shallow-curved  shells,  the  underlined  terms  are  ignored  and 
in  the  case  of  flat  plates,  the  radius  becomes  infinite  (R-*-®). 

Sanders'  complete  strain-displacement  expressions  are  valid 
when  (1)  the  Kirchof f-Love  hypothesis  holds,  (2)  middle  surface 
strains  and  rotations  out  of  the  middle  surface  are  small  in 
comparison  to  unity  (note:  the  displacement  components  are  not 
necessarily  small),  and  (3)  rotations  about  the  normals  to  the 
middle  surface  are  small  in  comparison  to  rotations  out  of  the 
middle  surface. 

To  allow  for  the  inclusion  of  transverse  strain  in  the 
variational  formulation.  Equations  1-3  are  used  to  define  the 
membrane  strains  in  each  individual  lamina,  and  the  lamina  summa¬ 
tions  process  provides  the  proper  strain-displacement  relations 
in  bending  (much  in  the  same  way  classical  lamination  theory  is 
developed).  This  approach  has  been  used  previously  in  Refer¬ 
ences  1,3~5  wherein  a  unique  model  for  inclusion  of  transverse 
shear  in  laminated  composites  has  been  shown  to  be  very  accurate 
in  the  prediction  of  initial  buckling  and  plate  bending  problems. 
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When  the  inplane  and  bending  strain-displacement  relations  are 
decoupled  in  this  manner,  the  Kirchof f -Love  hypothesis  is  effec¬ 
tively  removed  as  an  initial  assumption  and  the  transverse  shear 
strain-displacement  relations  become 


-  aw  *  uk+i  ~  uk 
Yzx  3x  t 

mk 


(*») 


V  —  v 

3w  k+1  k 

Yy  z  3y  t 

rak 


(5) 


2.3  Constitutive  Relations 

The  nonlinear  elastic  stress-strain  relations  used  (in 
material  coordinates)  are  based  on  those  of  Reference  1  and  are 
shown  below. 


(6) 


,°t  vtA 


+  nt^£t  + 

°At  U 


(7) 


At, 


=  b 


At 


At 


nAt,  A 
E  „ 


°t + 


K 


At  '■G 


t  T  At  >,n  At . 

( n  >  J  i 


At 


(8) 


^Gm 


m 


(9) 
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y 


zx 


k 


m 


k 


(10) 


Equations  9  and  10  reflect  the  transverse  shear  present  in  the 
model.  Note,  in  each  equation  a  Ramberg-Osgood  [6]  type  formula¬ 
tion  is  used  to  describe  the  inelastic  behavior  of  the  composite 
material.  These  relations  or  very  similar  ones  have  been  used  in 
the  solution  of  several  different  types  of  problems  [7,8],  and  in 
each  case,  the  resulting  numerical  model  has  been  shown  co  be 
quite  accurate  and  efficient. 


2.4  Modified  Hamilton/Reissner  Principle 

The  basis  of  the  analytical  procedure  developed  during  this 
research  is  the  extension  of  Hamilton's  principle  from  rigid  to 
deformable  bodies.  Hamilton's  principle  is  stated  as. 


<5 


(T  -  U  -  V)  dt 


0 


(11) 


where  t^  and  t^  are  arbitrary  time  values,  T  is  the  kinetic  energy 
and  U+V  is  the  total  potential  energy  of  the  system.  The  kinetic 
energy  is  the  sum  of  kinetic  energies  associated  with  each  of  the 
orthogonal  velocities  (axial,  tangential,  and  radial);  that  is. 


T  =  i. 
2 


a/2  b/2  h/2 

I  I  I 

-a/2  -b/2  -h/2 


+  +  (ffndzdydx  (12) 


3v  n  2 


3w  i  2 


where  p  is  the  mass  density  of  the  material.  Now,  if  Hamilton's 
principle  is  modified  by  replacing  the  strain  energy  U  with  the 
Reissner  functional  U"  [9],  the  variational  equation  becomes 
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6  -  /  (T  -  U")dt  =  0 

S 


(13) 


where  the  Reissner  functional  is 


a/2  b/2  h/2 

U"  =  /  /  /  P"  dzdydx 

-a/2  -b/2  -h/2 


(14) 


and  F"  Is  the  Reissner  functional 


F"  =  o  e  +  oe  +  t  y  +  t  y  +  t  Y  -  F '  „ 
xx  y  y  xy  xy  yz 'yz  zx'zx  fm 


(15) 


The  quantity  P'  is  the  stress-energy-density  function  which 
for  a  laminae  is  given  by 


p»  =  f '  +  f ' 

fm.  frn  rfrn 

k  ek  pk 


(16) 


where 


2  2 

>.  _  If0*  .  °t  (  V£t  .  Vt  Jl,  .  rn£,£t  .  n£ta, 

f™ek  '  ^  ^  '  (fT  (XT 


nt, It  +  nl t,t 


2 

Jit  ,  t  1  .  TJtt  ! 

WT*t  G  „  ,^k 


(17) 


=  | - i 

‘  (n  „+l 


\  (nrv,  s 


(n(*l)E(  ef  f 
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,  t  '  t  7  ,  2 

(nt+l)Et  eef  f 

.  lrl  1  1 

2  “  '  ’  (1+vti)i 

K»t  2 

<n*t+l,G*t  eff  lT“ 


(l-n  )/2 
1  it 


The  effective  strain  is  defined  by 


e 


eff 


k 


<v  + 


E»Et 


+  3 


2(1kt> 


2<1+nt> 


it 


it 


f  n  fc.  it+T]it.  i  N  .  , 

s;°v  ( 


nt. it+nit.t 


EtGn 


°*T  it 


H, 


1/2 


(18) 


For  transverse  shear 


(19) 
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where 


2  2 

{_Izz  +  I«i 

l2G  2G 

m  m 


and 


+  T 


zx 


(n  +l)/2 

)  m  }, 


(20) 


(21) 


where  Equations  6,  7,  and  8  apply  to  inplane  behavior  and  Equa¬ 
tions  9  and  10  are  applicable  to  transverse  shear  in  the  laminate. 
It  should  be  noted  that  the  derivative  of  the  quantity  F'  with 
respect  to  any  stress  component  provides  the  required  strain 
component. 

By  combining  the  strain-displacement  relations  and  consti¬ 
tutive  relations  of  Sections  2.2  and  2.3  with  the  modified 
Hamilton's  principle,  a  complete  new  variational  problem  is 
defined.  The  resulting  modified  Harailton/Reissner  functional  is 
presented  in  Appendix  A  in  vector-matrix  form  [10].  It  Is 
stated,  but  not  proven  herein,  that  the  first  variation  of  the 
modified  Hamllton/Reissner  functional  provides  both  the  governing 
equations  of  dynamic  equilibium  along  with  the  associated  stress- 
displacement  boundary  conditions  and  the  constitutive 
relationships . 


2.5  Temperature-Distribution  History 

For  each  kind  of  structure  examined  in  this  study,  the 
thermal  forcing  function  is  a  suddenly  applied  uniform  heat  flux 
acting  over  the  surface  z  =  -h/2  (see  Figure  1)  while  surface 
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z  =  h/2  is  presumed  to  be  perfectly  insulated.  These  conditions 
are  exactly  those  analyzed  by  Boley  and  Barber  [11],  Kraus  [12], 
and  Lu  and  Sun  [13]*  They  are  selected  here  for  convenient 
comparison  of  present  results  with  those  of  previous  analyses. 

The  temperature  distribution  history,  given  by  solution  of 
the  one-dimensional  Fourier  heat  conduction  equation,  is 


e(z,t)  = 


ail 

k 


r  <  1,1  z.2  1 

i,2t;  +  2  2  ~  h'  "6 


l  e-(n"'h)2t  cos  n,(i 


2 

it  n=  1 


n 


(22) 


where  q  is  the  heat  flux,  h  is  the  thickness  of  the  shell,  k  is 
the  thermal  conductivity,  and  k  is  the  thermal  diffusivity.  In 
Figure  2,  the  temperature-distribution  through  the  thickness 
corresponding  to  q  =  10  is  plotted  at  intervals  of  1/100  of  a 
thermal  period  for  a  duration  of  one  thermal  period. 


3.0  METHOD  OF  SOLUTION 

Extreraization  of  the  Reissner  functional  (Appendix  B) 
produces  the  governing  equilibrium  and  constitutive  equations  and 
the  consistent  natural  and/or  geometric  boundary  conditions. 

Once  these  governing  equations  are  obtained,  it  is  necessary  to 
effect  a  solution;  that  is,  determine  the  stresses  and  displace¬ 
ments  throughout  the  volumes.  The  exact  solution  of  the  nonlinear 
system  of  governing  equations  is,  in  general,  prohibitive.  Thus, 
approximation  methods  must  be  used.  This  being  the  case,  the 
extremization  of  the  modified  Hamllton/Reissner  functional  may  be 
used  in  more  direct  fashion  to  effect  an  approximate  solution. 
Since  the  vanishing  of  the  first  variation  of  the  modified 
Hamllton/Reissner  functonal  leads  to  the  governing  system  of 


10 


t  (SECONDS) 


Figure  2  Boley  and  Barber's  Histories  of  Thermoelastic 
and  Total  Moment  at  the  Center  of  the  Plate 
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equations  and  boundary  conditions,  then  it  is  possible  to  assume 
a  spatial  distribution  and  time  dependence  of  the  variational 
unknowns  satisfying  desired  boundary  conditions,  integrate  out 
the  spatial  dependence,  and  obtain  a  system  of  differential  equa¬ 
tions  for  the  time-dependent  variational  parameters  corresponding 
the  approximate  solution  of  the  governing  equations.  The  result¬ 
ing  system  of  nonlinear  time-dependent  differential  equations  may 
be  solved  by  any  number  of  numerical  procedures  (for  example, 
Newton-Raphson  for  converging  to  a  solution  at  any  particular 
time  step  and  Runge-Kutta  [14]  for  solution  in  the  time  domain). 
The  method  of  solution  is  outlined  in  a  subsequent  section. 

Application  of  the  Newton-Raphson  iterative  technique  results 
in  the  set  of  equations  contained  in  Appendix  C.  The  principal 
unknowns  are  the  stresses  (ii^, and  displacements 
(e,?,  and  g),  all  of  which  are  time  dependent.  Once  the  spatial 
variation  of  the  variational  unknowns  is  specified  (Section  3*1), 
the  Integrations  of  Appendices  B  and  C  can  be  performed.  The 
resulting  set  of  equations  are  syrabollically  written  as 


A*  =  X  t*  +  P  .  +  ft*  t* 

ext 


(23) 


where  the  individual  terras  are  expanded  in  Appendix  D. 


3.1  Displacement  and  Stress  Functions 

In  the  more  classical  procedures,  either  Minimum  Potential 
Energy  or  Complementary  Energy,  the  variational  unknowns  are 
either  the  displacements  or  the  stresses,  respectively.  The 
final  equations  involve  a  relatively  small  number  of  unknowns 
because  either  the  constitutive  relations  or  the  equilibrium 
equations  have  been  used  to  reduce  the  overall  number  of  governing 
equations  to  be  solved.  As  a  consequence,  the  final  equations, 
when  viewed  in  a  vector-matrix  form,  wherein  the  vector  of  unknown 
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quantities  are  either  displacement  or  stress  amplitudes,  have 
coupling  matrices  which  are  virtually  fully-populated  with  terms; 
i.e.,  the  coupling  is  strong.  However,  when  the  Reissner  varia¬ 
tional  approach  is  used,  the  vector  of  unknown  quantities  are 
both  displacements  and  stresses  representing  both  equilibrium  and 
stress-displacement  compatibility.  The  coupling  matrix  between 
the  variational  unknowns  is  observed  to  be  very  sparse;  i.e., 
loosely  coupled. 

This  loose  coupling  is  a  problem  because  an  inappropriate 
choice  of  displacement  and  stress  functions  can  lead  to  a  situa¬ 
tion  wherein  the  computed  integral  that  couples  the  physical 
terms  can  be  zero — this  would  imply  that  the  computational  proce¬ 
dure  would  not  recognize  the  physical  coupling  that  is  present  in 
the  real-life  structure.  As  a  consequence  of  this  mathematical 
consideration,  it  is  necessary  to  select  a  complete  set  of  func¬ 
tions  (displacement  and  stress)  that  admit  coupling  between 
physical  quantities  that  should  be  coupled. 

It  should  also  be  noted  that  it  is  possible  to  select 
functions  that  provide  for  coupling  when  it  should  not  occur. 

This  can  also  lead  to  fictitious  results.  In  any  event,  the 
first  and  foremost  requirement  to  be  placed  on  the  selected 
displacement  and  stress  functions  is  that  the  proper  kind  of 
coupling  be  maintained  between  the  physical  entities  of 
displacement  and  stress. 

The  next  requirement  to  be  satisfied  by  the  distribution 
functions  is  that  they  represent  the  physical  behavior  of  the 
structure  under  analysis.  In  practice  this  is  usually  not  too 
difficult  in  the  sense  that  It  is  usually  easy  to  select  either 
trigonometric  or  polynomlc  functions  that  match  approximately  the 
expected  displacements  and  stresses  of  the  structure  to  be 
analyzed.  In  this  sense,  for  complicated  problems,  wherein  the 
exact  form  of  the  distributions  is  not  known  (as  with  the  present 
problem),  it  is  necessary  to  select  several  "generic"  forms  and 
experiment  with  the  combined  distributions  to  determine  which 
ones  provide  viable  solutions  without  disturbing  the  need  for 
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proper  coupling  as  discussed  previously.  At  this  point,  the 
selection  and  verification  of  displacement  and  stress  functions 
used  to  solve  the  problem  of  the  buckling,  postbuckling,  and 
crippling  of  laminated  composite  shallow  shells  under  combined 
axial  compression  and  shear  [5]  took  approximately  seven  months 
of  intense  effort  in  spite  of  the  fact  that  the  initial  selection 
of  trial  functions  was  achieved  in  approximately  two  weeks.  The 
issue  of  loose  coupling  was  a  major  factor  in  the  time  required. 

Although  not  mentioned  yet,  but  just  as  important,  there  are 
additional  requirements  that  are  placed  on  both  displacement  and 
stress  functions.  For  instance,  displacement  functions  must 
satisfy  the  boundary  conditions  of  the  problem  under  consideration. 
Stress  functions  must  satisfy  equilibrium  across  cross-sections 
as  a  minimum  (based  upon  the  experience  of  the  principal  investi¬ 
gator)  ,  although  this  requirement  can  be  waived  under  the  proper 
circumstances . 

It  should  be  remembered  that  in  using  the  Relssner  varia¬ 
tional  principle,  the  constitutive  relations  are  not  satisfied 
concurrently  at  every  point  in  the  structure.  At  first  thought, 
this  rather  pointed  observation  is  disturbing  because  most  knowl¬ 
edgeable  enlgneers  familiar  with  mechanics  would  state  that  the 
constitutive  relation  must  be  satisfied  at  every  point  for  the 
solution  to  be  considered  correct.  In  actual  practice  this 
presents  no  particular  problem  if  some  ordering  of  priorities  is 
kept  In  mind. 

First,  the  objective  in  using  the  Reissner  variational 
approach  is  to  obtain  an  accurate  knowledge  of  the  state  of 
displacement  in  the  structure  at  the  minimal  computational  cost. 
This  Is  usually  achieved  by  careful  selection  of  a  complete  but 
modest  number  of  displacement  functions  which  satisfy  the  boundary 
conditions.  Stress  functions  are  then  cnosen  to  satisfy  equi¬ 
librium  across  some  cross-section  or  area,  but  not  at  every 
point.  These  same  stress  functions  must  couple  correctly  with 
the  displacement  functions  in  such  a  way  that  the  finally  com¬ 
puted  displacement  amplitudes  and  distributions  are  accurate.  In 
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actual  practice  It  takes  a  combination  of  knowledge,  intuition, 
and  perseverance  to  achieve  this  end--but  it  is  possible,  as 
documented  in  earlier  similar  efforts  [3,4,5]. 

In  the  end,  when  the  numerical  results  are  produced,  the 
stress  amplitudes  are  essentially  discarded  in  the  sense  that  the 
numerical  values  of  stress  may  or  may  not  be  accurate--it  is  not 
important  because  the  goal  of  the  approach  is  to  obtain  accurate 
information  about  displacement.  Prom  displacements  the  strains 
may  be  readily  calculated  because  the  strain-displacement  rela¬ 
tions  are  in  general  accurate  and  a  single  derivative  does  not 
lead  to  gross  inaccuracies.  Stresses,  of  an  accurate  nature,  can 
be  obtained  via  the  use  of  the  original  complete  set  of  nonlinear 
constitutive  relations. 

The  discussions  of  the  preceding  paragraphs  give  an  idea  of 
the  difficulties  associated  with  the  selection  of  displacement 
and  stress  functions  for  use  in  the  Reissner  variational  principle. 
In  general,  these  comments  apply  to  any  particular  problem;  how¬ 
ever,  for  this  particular  research  concerned  with  the  effects  of 
delaminations,  there  are  some  additional  considerations  that  must 
be  accounted  for.  In  addition  to  the  unique  application  of  the 
Reissner  and  Hamilton  principles,  the  use  of  the  direct  varia¬ 
tional  approach  to  solve  a  problem  which  inherently  couples  local 
and  global  displacement/stress  fields  Is  also  noted  as  a  signifi¬ 
cant  challenge. 

It  is  anticipated  that  interactions  between  local  and  global 
displacement/stress  fields  can  occur  deleteriously ,  as  discussed 
previously.  Thus,  selection  of  global  and  local  functions  is 
potentially  a  difficult  and  time  consuming  process  even  though 
the  starting  functions  seem  intuitively  obvious;  i.e.,  the  local 
stress  fields  near  a  crack-tip  representing  a  delaraination  are 
generally  known. 
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3.1.1  Contiguous  Composite 

Based  on  the  earlier  work  of  Stroud  and  Mayers  [15] >  the 
displacements  are  taken  as 


u  =  u1 c  (24) 

v  =  v^n  +  v 2 ( 1  -  C2)(n  -  n^)  (25) 

w  =  W22(1“C2)  (1-n2)  +  w2^(l-C2)  (  n2-^)  +  wl42  (  C2-^)  (l-^2) 

(26) 


where  c  and  n  are  dimensionless  goeraetric  coordinate  parameters 
defined  as 

C  =  2x/a  (27) 


and 


n  =  2y/b  (28) 

For  flat  plates,  the  second  term  for  the  expression  for  v  is 
omitted . 

It  is  apparent  that  for  the  linear  thermoelastic  response  of 
a  plate  with  unrestrained  edges,  the  t>  tions  assumed  for  inplane 
displacements  are  exact.  The  transverse-displacement  expression 
is  a  fourth-order  polynomial  with  all  except  three  of  the  coeffi¬ 
cients  eliminated  by  enforcing  geometric  boundary  conditions  and 
symmetry  requirements  (that  is,  vanishing  displacement  at  the 
supports  and  symmetry  about  the  coordinate  axes). 

The  distributions  of  bending  and  twisting  moments  assumed 
for  the  thermoelastic  solution  are 


M 


xB 


n  r  i  o  p 

Nt  +  [  (1-v)  (Mt  -  —  Nt)  -  Mx22](1-c  )n 
H  H 
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+  Vix22^1  ~  C2(1-n10)]  +  Mx21<  (l-;2)  ( n2-ni<)  +  Mx2|2  ( c2-;11 )  (1-n2) 

(29) 

7^  7*  p  1  f) 

MyB  =  "MT  +  NT  +  [(1'^(Mt  "  Nt)  -  My22]C  (1-n  ) 

+  My22[l-(1-C10)n2]  +  My24(^2-;4)(l-n2)  +  My42(l-i;2)(n2-r14) 

(30) 

MxyB  =  Mxyll  3ln  i  «  3ln  f  "  *  Mxyl3  Sln  f  t  sln  T-  " 

*  Wjyj!  sln  2^  c  sln  2  "  (3D 


3.1.2  Delaminated  Composite 

For  a  delaminated  composite  specimen  the  out-of-plane  dis¬ 
placement  w  Is,  In  general,  a  function  of  the  through  thickness 
coordinate,  z,  as  well  as  the  other  spatial  coordinates  and  time; 
that  Is,  w  =  w(z,x,y,t).  Equivalently,  the  out-of-plane  displace¬ 
ment  is  considered  to  be  a  function  of  the  lamina  identification, 
k;  thus,  each  lamina  can  have  a  different  functional  form  for  the 
out-of-plane  displacement.  For  Instance,  in  the  local  region 
around  the  delaminatlon,  a  potential  choice  for  displacement 
function  Is 
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“local  =  “dl1  +  003 

where  w^  is  the  amplitude  of  the  delaraination,  x  and  y  are 
spatial  coordinates  located  at  the  center  of  the  assumed  circular 
delamination  of  radius  a,  and 

w  =  w  .  ,  _  +  w..  _ 

global  local 


where  w 


global 


was , 


Corresponding 


for  a  contiguous  composite,  the  same  as  w. 
to  this  out-of-plane  displacement  are  local 


stress  fields  analogous  to  those  near  a  crack  tip  (see  Sih  [16], 


for  instance) . 


3.2  Iterative  Solution  Technique  in  the  Time  Domain 

When,  at  any  one  time  step  subsequent  to  updating  material 
properties  and  temperatures,  the  governing  equations  are  solved, 
then  +  0  and  the  resulting  differential  equation  is 

®t*  +  +  P  =  0  (32) 

GX  u 

Defining 

y  =  i*  (33) 

then 

+  ffit  +  P  «_  *  0  (3*0 

ext 
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Equations  (33)  and  (3*0  can  be  written  as  a  set  of  first  order 
differential  equations: 


■— 

" m 

( .  \ 

““ 

/  \ 

*"  “ 

flf 

0 

< 

-*• 

y 

A 

0 

< 

+ 

X 

>  + 

^ext 

= 

3 
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\  J 

0  - 

r 

J, 

0 

3 

The  above  equation  can  be  solved  with  the  use  of  the  Runge- 
Kutta  [14]  approach. 


4.0  RESULTS  AND  CONCLUSIONS 

The  objective  of  this  Phase  I  research  was  to  derive  the 
equations  governing  the  response  of  a  delaminated  composite  to  a 
rapid  thermal  pulse.  The  approach  taken  wa3  one  of  developing  a 
variational  formulation  based  on  a  unique  coupling  of  Hamilton’s 
and  Reissner's  variational  theorems.  The  resulting  equations 
have  been  cast  into  a  form  that  is  suitable  for  programming  on  a 
digital  computer.  The  solution  of  these  equations  provides  a 
description  of  the  state  of  displacement  and  stress  within  the 
delaminated  composite. 
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APPENDIX  A 


MODIFIED-HAMILTON/REISSNER  FUNCTIONAL 
IN  VECTOR-MATRIX  NOTATION 
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Hamilton's  principle  is  stated  as 


t2  t2 

6  /  (T  -  U  -  V)dt  =6  /  (T  -  U" )dt  =  0 

t  t 

1  1 


(A-l) 


where  the  Reissner  functional  in  vector-matrix  notation  [10]  i3. 


U"  =  f  f  x^  p^  t  {f  —  k  +  +  Q  -2-S  e} 

x  Jy  1  1  A1  rl  rfmu  2  s  9x  3x  g  6±  wx  3x  e) 


+  ixT  XT  PT  t  ft  ±  2T  11  ?  +  3  li  ? 

n2  2  2  fm1  2  fi  3y  3y  g  3y 


~T  «T  ~T  f  1  *T  ~T  xT  x  «  ~  £  3  y  *  1 

-  h2  X2  R2  g  tfm  Y  R-  +  f  e  Qy  tfm  P2  X2  h2  g  R 


*T  ~T  ~T  ~  ~  -*■  *T  ~T  ~T  1 
+  31  Qt  P,  X0  h0  P  e1  Qi  Ar} 


y  2  2  2 


y  R2 


+  hT  xT  PT  t  ft  gT(^  H  £)  +  Q  —  e  +  Q  ? 
n3  3  3  fmu  6  l3x  3y  g;  wx  3y  e  wy  3x 


*  ^  ll  vtT  rT  dT  P  n  a 

g  3X  h3  x3  P3  fcfm  Qy  0  *} 


+  aT  f  rJ{!  g  +  B  Q~x  S  e} 


+  ST  rJ{I  |^-  2  +  B  5y  s  f} 


3  3 


=  1  M»1  ^  ^LM  AL  hL  ~  VTWl  ”  ayAT?2A2h2 


24 


Anamet  Laboratories,  Inc. 


HAYWARD,  CALIFORNIA 


VT  ~T  g-T  ~ 

-a  e  R  G  R  e  a 

x  in  x 


4-T  ~  ~  -* 

b  n  R  Q  R  n  b 
y  m  y 


3  3  +T  ~T  ~T  N  3  3  * 

l  l  hj  Aw  Pw  l  l  1  f1T  ( 2  6,  ,  ,  ) 

L=1  Mil  M  M  M  k=l  1=1  3=1  iLk  1J'1 


3  3 

_  "d  ->fp  mJV  rs/rP  />.•  /s,  *s/  =*>  _  l/-  ^  > 

[  [  I  h.  a;  p:  d  .  ,  .  p  x  hi  K  f,  pT  AThT 
£=1  m=l  £  4  A  £miJk  m  m  m  jmk  L  L  L 


(N-l) 


y  UT  £T  Bk  Rx£  a  +  ST  nT  fp  Bk  Ry  n  S)  4}dxdy 
1 


y  "k  y 


( A-2) 


and  similarly,  the  kinetic  energy  is 
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FIRST  VARIATION  OF  THE  MODIFIED-HAMILTON/REISSNER  FUNCTIONAL 
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The  governing  equations  which  result  from  taking  the  first 
variation  of  the  modified  Hamilton/Reissner  functional  (Eq.  (A-l)) 
are  as  follows: 
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NEWTON-RAPHSON  ITERATIVE  SOLUTION  TECHNIQUE 
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Application  of  the  Newton-Raphson  iterative  technique  provides 
the  following  equations: 
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